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aircraft control vector 


A concept for automatic terminal-area guidance, 
comprising two modes of operation, has been devel- 
oped and evaluated in flight tests. In the first 
or predictive mode, fuel-efficient approach trajec- 
tories are synthesized in fast time. In the second 
or tracking mode, the synthesized trajectories are 
reconstructed and tracked automatically. An energy 
rate performance model derived from the lift, drag, 
and propulsion-system characteristics of the air- 
craft is used in the synthesis algorithm. The 
method optimizes the trajectory fo*- the initial air- 
craft position and wind and temperature profiles 
encountered during each landing approach. The 
paper describes the design theory and discusses the 
results of simulations and flight tests using the 
Augmentor Wing Jet STOL Research Aircraft. 
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* airspeed, ft/sec or knots 

- final and initial airspeeds of air- 
craft, respectively, ft/sec or knots 

- wind speed in direction of ground 
track, knots 

- aircraft weight, lb 

■ perturbation state vector 

* final and initial x coordinates 
of aircraft, respectively, ft 

■ final and initial y coordinates 
of aircraft, respectively, ft 


List of Symbols 


angle of attack, deg 


drag force, lb 


■ inertial flight-path angle, deg 


V d f 


distance of backward and forward 
integration, respectively 


* Aerodynamic flight-path angle, 
rad or deg 


cruise distance, ft 

length of ground track from initial 
to final position of aircraft, ft 


airspeed rate correction due to 
wind shear, ft/sec 2 

crosstrack error, ft 


■ energy, ft 


crosstrack error rate, ft/sec 


energy rate, ft /sec 


flap angle, deg 


E E 

nmax* nmin 


maximum and minimum available 
energy rate, respectively, ft/sec 

perturbation state and control 
distribution matrices, respectively 

acceleration of gravity, ft/sec 2 

altitude, ft 


maximum flap angle, deg 

fraction of energy rate used for 
changing speed 

command pitch angle, deg 

vectored thrust, in degrees of 
nozzle angle 


V h i 


It k ♦ 

♦y f %y 


• final and initial altitudes of air- 
craft, respectively, ft 

■ feedback gain matrix 

• lateral error and error rate feed- 
back gains 

• lift force, lb 


v *i 


« throttle setting, in percent RFM 

■ fraction of available energy rate 

■ commanded and reference bank angles, 
respectively, deg 

- final and initial ground heading of 
aircraft, respectively, deg 


S * speed along ground track, ft /sec 

T * thrust force, lb 

t • time, sec 

u * perturbation control vector 

*Research Scientist. Member A1AA. 

This paper is declared a work of the U.S. Government and 
therefore is in the public domain. 


I ntroduction 

In the past, terminal-area guidance system 
design for aircraft has concentrated primarily on 
automatic glide slope tracking, flare, and touch- 
down. During recent years, designs have been devel- 
oped to provide automatic guidance along curved and 
decelerating approach paths. * This increased capa- 
bility was made possible through the integration of 
digital computers into the flight guidance system. 
However, even In the more advanced designs, automatic 
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guidance is limited to 3 few prestored three- 
dimensional flight paths, as in Ref. 1. While the 
ability to fly complex prestored trajectories is 
essential, it cannot give optimum performance under 
actual terminal-area operating conditions as shall 
be explained. 

First, a prestored trajectory cannot optimize 
fuel consumption or a similar performance measure 
under actual operating conditions. Optimum trajec- 
tories depend significantly on aircraft gross 
weight, wind and temperature profiles, and on the 
initial state of the aircraft. These variables 
cannot be predicted with the required precision 
prior to takeoff. To prestore optimum trajectories 
for each of the conditions likely to be encountered 
would result in an impossibly large memory require- 
ment. Therefore, prestored trajectories must neces- 
sarily represent a compromise in performance. 

Second, in existing systems the pilot must fly 
the aircraft manually from its current position to 
the starting point of the trajectory. This flight 
segment is known as the capturing maneuver. Three- 
dimensional, curved trajectories can be difficult to 
capture manually, and, if the trajectory alto 
includes a specification of landing time, as is the 
case in 4D guidance, the capturing maneuver cannot 
be done by the pilot without computer assistance. 
Therefore, the capturing maneuver, because of its 
variability, can only be generated by onboard 
trajectory synthesis. 

Third, aircraft in high density airspace are 
usually controlled by air traffic control vectors 
and during this period cannot follow a prestored 
flight path. Synthesis of a trajectory can only 
begin after the aircraft has received its final 
vector and has been cleared for approach. But the 
initial position of the aircraft at that time varies 
between approaches, thus trajections require 
onboard synthesis. 

An initial design of a four-dimensional guidance 
system embodying the concept of onboard trajectory 
synthesis, including an advanced capture law, was 
aeveloped and flight tested onboard a Convair 340 
aircraft equipped with the STOLAND avionics. 2 In 
the design described here, horizontal trajectories 
are generated by the method of Reference 2, but 
vertical and speed profiles are synthesized using 
simplified aero/propulsion performance models of 
the aircraft. This results in profiles that are 
optimum for fuel conservation. Design of the con- 
trol law for tracking the synthesized trajectory is 
based on the linearized perturbation guidance 
approach. Since the perturbation equations are air- 
craft configuration-dependent, gain scheduling is 
used in the feedback law. 

The Augmentor Wing Jet STOl Research Aircraft 
(AWJSRA) was chosen as the test vehicle for this 
concept. This type of pcwered-llft aircraft is 
highly cost-sensitive to operational procedures in 
the terminal area. It also exemplifies particularly 
well the unique problems of powered-lift aircraft, 
namely, high fuel consumption in the STOL mode; 
dependence of both lift and drag on thrust; and an 
excess of controls over the minimum number needed to 
determine path and speed. These factors suggest 
that trajectory optimization could greatly increase 
the operational efficiency of the aircraft. Imple- 
mentation of this concept was facilitated by the 


existing installation of the STOLAND avionics system 
onboard the aircraft. 

Energy Rate Model and Selection of R c ferem. e_ Con t nls 

An energy rate model of aircraft performance 
has been found to yield a compact and sufficiently 
accurate representation of performance for terminal- 
area trajectory synthesis. In this section a per- 
formance model based on energy rate is derived and 
then applied to determine the optimum reference 
controls* for synthesizing trajectories. 

Consider the standard expression for energy 
rate written as 


where 


dE 

dt 


(T - D)V a 
W 


( 1 ) 


E * h + \ 2 ( 2 ) 

with constraint L * W (Ref. 3). It is assumed 
throughout this paper that flight-path angles are 
small such that cos y a 1 and sin y a * y a . Fur- 
thermore, it is assumed that flight-path angle rates 
are so small that their effect on lift is negligible. 
Differentiation of Eq. (2) with respect to time 
gives an equivalent expression for energy rate: 


dE 

dt 


dh + i Va _a 

dt g d dt 


(3) 


Equations (1) and (3) can be nondimensionalized by 
dividing them both by V a . The resulting quantity 
on the left side, (1/V a ) (dE/dt) , is defined as the 
normalized energy r-te £ n , or energy rate for 
short. By using the relation (dh/dt) % v 3 Y a * 
the two relations for £ n become 


E 

n 



1 dV a 

Y a + IdT 


with constraint L ■ W. 


(4) 

(5) 


Equation (4) specifies the energy r^te as a 
function of the difference between thrust and drag, 
subject to the constraint that lift equal weight. 
Thrust and drag are in turn functions of the controls 
producing forces in the flight-path direcUon, namely 
throttle tt, flap angle 5f, nozzle angle v 
(vectored thrust), and angle of attack a. Equa- 
tion (5) determines the relationship between flight- 
path angle and acceleration for the energy rate 
calculated from Eq. (4). Equation (5) indicates 
that, in particular, a giver, energy rate may be 
utilized to fly at flight-path angle y a with con- 
stant airspeed, or to fly at zero flight-path angle 
with acceleration rV a /dt. An infinity of other 
combinations of y 9 and dV a /dt can also be chosen 
to yield the same energy rate. This makes possible 
a simplifying dichotomy in the trajectory synthesis, 
namely, at any time the desired energy rate Is 
selected first by choice of appropriate controls and 
then the linearly related quantities of y a and 
dV a /dt are selected to generate the specifics of 
the flight path. The next section develops the 
complete synthesis algorithm based on this approach. 
Here we elaborate on the determination of the 
functional dependence of energy rate on the force- 
producing controls. 


\ 
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Since the STOL aircraft studied in this paper 
has four controls to achieve a specified energy rate 
and to maintain lift equal to weight, there is an 
excess of two controls over the minimum number 
needed for a simultaneous solution to Eq. (A) and 
the constraint L ■ W. These two extra degrees of 
freedom in the controls are exploited to minimize 
power setting and, therefore, fuel flow at every 
energy rate. This optimization problem is restated 
in equivalent form as the maximization of energy 
rate for a given power setting: 


E n 00 - 

T - D 

max u 

v,a,6 f 

(6) 

Constraint : 

: L(tt , v ,a , Sf ) - W 

(7) 


The maximization must obey various inequality con- 
straints on the controls: 

-10.5° < a < 19.5° 

6° < v < 100° 

5 ’ 6 ° < 6 f * fmax (V a ) < Fla P P^ard] 

In addition, a lift or maneuver margin must be 
satisfied at every point to guarantee sufficient 
normal force for changing the flight path. Pilots 
familiar with this aircraft specify that at least 
0.4 g of normal acceleration must be attainable at 
any time by an increase in the angle of attack 
alone. 

The use of Eq. (6) results in the selection of 
the controls that yield the maximum attainable 
energy rate at each thrust setting. This ensures 
the efficient use of thrust at any energy rate that 
requires more than the minimum thrust. But energy 
rates more negative than those attainable Dy Eq. (6) 
are also of interest. Such negative er.ergy rates 
must occur at the greater of minimum or idle thrusts 
required by the maneuver margin. At a particular 
airspeed, a decrease in the energy rate below the 
minimum attained througn Eq. (6) can be effected by 
increasing the vectored thrust angle v and/or the 
flap angle <5f. The third control, angle of attack 
a, is needed to satisfy thy constraint L * W. The 
two degrees of freedom in the controls can be 
exploited to minimize noise exposure along the 
ground crack. Noise under the aircraft is known to 
increase as the nozzles producing the vectored 
thrust are turned downward. Therefore, a further 
decrease in energy rate is achieved by first increas- 
ing flap angle until it reaches its limit or placard 
value and only then by Increasing nozzle angle. 

The result of applying these procedures to the 
AWJSRA is shown in Fig. 1 for a weight of 38,000 lb, 
sea-level altitude, and standard temperature. The 
figure gives the envelope of energy rate vs indicated 
airspeed with throttle, flaps, and vectoring nozzle 
as parameters. Angle of attack Is not plotted to 
avoid cluttering the figure. At any airspeed, the 
^nmax an< * Enmin curves define the range of per- 
missible energy rates. The optimum controls for a 
given airspeed and energy rate are determined by 
Interpolation between contours of constant controls. 
For example, at an airspeed of 105 knots and 
L. * -0.17, the optimum controls are found to be: 

Gf - 26*, v - 6*, ir - 84% (point A, Fig. 1). Angle 
of attack (not ahowr) is 8.4*. Maximum energy rate 
with minimum thrust occurs at 112 knots (point B) 
and corresponds approximately to (L/D) max « 10* 


It should be noted that the force-producing 
controls in this experimental STOL aircraft have 
unusual characteristics that account for the rela- 
tive complexity of Fig. 1. Throttle affects both 
lift and drag at all speeds, but the effect on lift 
is greatest in the STOL regime below about 80 knots. 
The thrust magnitude produced by the vectoring 
nozzle, referred to as the hot thrust, is also 
coi * rolled by the throttle and accounts tor about 
60% of the total thrust produced by the two engines. 
The remaining 40% of the thrust which is the cold 
thrust produced by the fans, energizes the augraentor 
wing to increase lift at STOL speeds. 

The relationship between the controls and the 
energy rate is revealed more clearly in Fig. (2) 
at the example airspeed of 105 knots. Many such 
plots at various airspeeds would be required to 
illustrate the complete dependence of the controls 
on energy rate. As the energy rate decreases below 
Its maximum value of 0.28, throttle decreases nearly 
linearly until idle throttle is reached. In this 
interval flaps Increase only slightly while nozzle 
angle remains at minimum and angle of attack 
increases. At more negative energy rates, flaps 
become the dominant control until they reach the 
placard value of 40° at this airspeed. Angle of 
attack decreases sharply as flap angle increases. 
Finally, nozzle angle increases toward its maximum 
value 3f 100® as the energy rate decreases toward 
its negative limit of -0.3. 

In the flight implementation of the algorithm, 
four diagrams as shown in Fig. 1 are utilized, two 
for sea-level altitude at weights of 38,000 and 
48,000 lb and two others for 5000- ft altitude at 
similar weights. Experience indicates that these 
are sufficient data to interpolate the controls 
adequately. Each diagram requires 124 words of 
memory in the airborne computer. The small circles 
in Fig. 1 indicate the locations of points that are 
stored. The energy rate data are also corrected 
for deviations from the standard temperature profile. 
Correction is done by computing a thrust setting 
corrected for temperature deviations. 

Synthesis of Complete Profiles 

In the preceding section rhe criteria of fuel 
conservation and noise reduction were used to deter- 
mine the four refei ..ce controls of throttle, noz2le 
angle, flap angle, and angle of attack as functions 
of the energy rate. This approach replaced the 
problem of selecting four control variables with 
the simpler problem of selecting a single, equivalent 
variable, namely, the energy rate. In this section 
we make use of the energy rate variable in generat- 
ing efficient terminal-area trajectories. 

Tha problem of terminal-area-trajectory syn- 
thesis can be stated as the specification of rules 
for flying an aircraft with initial state vector 
I x i> y±» ^i» v ai3 to a final state vector 

[x f , y f , hf, V a f]. To be of practical interest, 
such rules must generate efficient and flyable 
trajectories connecting various initial and final 
state vectors. By specifying a performance criterion 
such as fuel consumption, we can fit this problem 
into the framework of optimal control theory. 

However, the difficulty of solving an optimal control 
problem characterized by a five-element state vector 
makes this approach computationally Impractical for 
in-flight implementation. Following Ref. 4, we have 
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adopted the simplifying procedure of separating the 
synthesis problem into two essentially independent 
problems . 

The first problem consists of synthesizing the 
horizontal or 2D trajectory. References 4 and 5 
give several algorithms for computing near-minimum- 
distance 2D trajectories as s sequence of an initial 
constant radius turn, straight flight, and a final 
constant radius turn, where the turn radii are 
chosen so as to avoid exceeding a specified maximum 
bank angle at the maximum ground speed encountered 
in each turn. A description of the algorithm used 
in the flight implementation can be found in 
Ref. (5). Figure 3 illustrates the 2D trajectories 
computed by the algorithm for several initial posi- 
tions, Pi, in the terminal area. Note that the 
terminal point, Pf, lies on an extension of the 
runway centerline, and that the heading angle 
of all trajectories is equal to the runway heading 
at that point. Thus, the algorithm always generates 
2D trajectories that match the initial and final 
state vector components x*, y*, and x f , y^, . 

The second problem, solved after the horizontal 
trajectory has been computed, consists of cynthesiz- 
ing efficient speed and altitude profiles which 
match the initial and final speeds and altitudes 
v i» hi* and Vf, hf, respectively. 

The horizontal distance of the trajectory dh, 
a known quantity computed in the previous step, adds 
a third boundary condition to be satisfied by the 
profiles. While this three-state optimal-control 
problem is much simpler to solve than the original 
five-state problem, it is still too complex for 
onboard-computer implementation. A simpler algo- 
rithm was therefore developed that generates near- 
Qneed-al fitude profiles by matching the 
general characteristics of optimum fuel and noise 
trajectories studied in Refs. 6 and 7, respectively. 
We briefly explain the rationale for this algorithm 
with reference to descent, which is the most dif- 
ficult case. 

It was found in Ref. 6 that the descent portion 
of a minimum fuel descent trajectory is characterized 
by a delay in the start of the energy decrease as 
long as possible consistent with meeting end con- 
straints of speed and altitude. Furthermore, the 
energy change consists initially of descent to the 
final altitude at near-constant indicated airspeed 
followed by a rapid airspeed deceleration In level 
flight. Most of the energy change takes place at 
minimum throttle, as one might expect for minimum 
fuel flight* Minimum noise descent, profiles com- 
puted in Ref. 7 are similar in that they also delay 
the start of energy decrease as long as possible, 
but they approach the final altitude in a steep 
descent to maximize the aircraft's altitude above 
the ground near the runway. This means that the 
deceleration to the final airspeed takes place 
before the start of descent or during the early 
portion of the descent* Thus the two types of 
descent profiles differ prirarily in the way they 
proportion the use of available energy rate to 
decrease altitude and airspeed. 

To facilitate the synthesis of such profiles, a 
family of decreasing (and by extension, increasing) 
energy profiles, which Include the two types 
described as special cases, is defined by two param- 
eters, o and c. The first parameter, o, selects the 
fraction cf minimum /maximum available energy rate. 


Erunin> (^runax) t0 be use d f° r decreasing/ increasing 
energy. The values of E nmin and E nmax can be read 
from Fig. 1 af each indicated airspeed. The second 
parameter, t, determines the fraction of the selected 
energy rate to be used for decelerat ion/acoelerat u»n. 
Then, for particular cnoices of * and c, the energy 
rate, airspeed, flight-path angle, altitude, and 
horizontal distance are computed as follows: 


E 

= oE 

0 

s':**! 

(7) 

n 

nmin 




V 

a 

= gcE n 

0 

< c < 1 

(8) 

y 

» (1 - e)E 



(9) 

a 

n 




. 

* V y 



(10) 

n 

a a 




S 

= V cos y 
a a 

-t- 

V 

w 

(11) 


where V w is the along-track component of wind 
speed. Note that Eqs. (7)-(9) are consistent with 
Eqs. (4) and (5) for all values of o and e. 
Decreasing/ increasing energy profiles are generated 
by integrating Eqs. (8), (10), and (11) for parti- 
cular choices of o and c. 

To illustrate the effect of the parameter e 
on the descent /decelerat ion profiles, assume 

» -0.13, independent of speed, and let the air- 
speed to be achieved at touchdown be 100 ft /sec. To 
achieve the desired boundary conditions, Eqs. (8), 
(10), and (11) are integrated in backward time 
starting with the speed and altitude at touchdown. 

The resulting airspeed and altitude profiles are 
plotted as a function of distance to touchdown in 
Fig. 4 for e - 1, 0.5, 0.0. The profile for e - 1 
is seen to approximate the minimum fuel, for e - 0, 
the minimum noise descent, and for e = 0.5, a 
compromise between fuel and noise minimization. 

To achieve minimum fuel or noise performance, 
changes in energy should be made at maximum rate. 

This is accomplished by setting o to unity and 
thereby following the E nm in contour during descent 
and deceleration. However, for the aircraft under 
study this yields energy rates too negative for safe 
operation in the terminal area at some airspeeds. 

A limit less than one is also necessary to reserve 
energy rate for perturbation control. A practical 
upper limit on o is about 0.9 for the AWJSRA. In 
the flight implementation, the two profile parameters 
are keyboard entries that allow the pilot to choose 
values appropriate for each landing approach. In 
addition, the pilot can specify the maximum decelera- 
tion and descent angles via keyboard entry. The 
maximum safe deceleration for this aircraft is 
limited to about 0.06 g by the maximum rate at 
which flaps can be extended. The synthesis algorithm 
is configured to decrease a below its limit if that 
is necessary to satisfy these constraints. 

The backward time integration described above 
generates an increasing (in backward time) energy- 
profile starting at the desired final speed and 
altitude. To complete the synthesis of the descent 
trajectory we still need rules for matching this 
profile to the initial speed and altitude of the 
aircraft* The freedom of the aircraft to maneuver 
in altitude is restricted by air traffic control as 
well as passenger comfort considerations. Thus, as 
an aircraft approaches a terminal area. It Is 
generally not allowed to climb above its initial 
approach altitude for the purpose of optimizing the 
approach trajectory* The aircraft must hold this 
altitude until starting the final descent. However, 
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while flying at altitude hi, it may change to a 
new airspeed, V at called the terminal-area speed 
which can be higher or lower than the initial speed 
V a i- Unless specified by the pilot via keyboard 
entry, it is chosen to minimize fuel use per unit 
distance, and is 140 knocs for this aircraft (it 
would be 220-250 knots for conventional jet 
transports) . 

The various rules contained in the preceding 
two paragraphs can now be combined to yield the 
complete algorithm. The synthesis begins with the 
backward time integration from final conditions 
hf, v af using the specified o and e. If the 
altitude reaches its target value of before the 

airspeed reaches its target value of V at , we set 
e ■ 1 and then continue Lhe backward time integra- 
tion until the airspeed has also achieved its target 
value. When setting e * 1, the flight-path angle 
is forced to zero and the energy rate is used 
entirely for accelerating (in backward time) toward 
V at . On the other hand, if the airspeed reaches its 
target value before the altitude does, we set e * 0. 
This stops the airspeed change and uses the energy 
rate entirely for increasing the altitude toward its 
target value of h^. When the second and last vari- 
able reaches its target value, we set a - 0, i.e., 

En * 0, thus completing the backward time integra- 
tion. Next, we begin a forward time integration to 
get the distance required to change speed from V a ^ 
to V at with e * 1. Let the distances for the 
backward and forward Integrations be d^ and df, 
respectively. A valid trajectory has been generated 
if the cruise distance d c , computed from 


is nonnegative, i.e., d c > 0. If d c is negative, 
the synthesis has failed because the aircraft is too 
close to the capture point Pf. 

Figure 5 illustrates the various segments of an 
approach trajectory synthesized by the algorithm. 

As before, we assume for simplicity that 
E|| - -0.13, a constant. Other parameters defining 
the problem are Indicated in th«* figure. Note that 
the initial descent at y a * -7.5 shallows to 
y a • -3.75 to allow the aircraft to decelerate. The 
reference controls for this trajectory can be inter- 
polated from Fig. 1. 

We conclude this section by mentioning briefly 
other important features of the algorithm included 
in the flight implementation. The airspeed decelera- 
tion is corrected for known wind shears, which are 
computed frem a knowledge of V w (h) if available. 

The wind shear correction factor is 

AV a - -(dV w /dh)V a y a (13) 

and is added to the right side of Eq. (8) to obtain 
the corrected airspeed rate. Furthermore, the refer- 
ence controle are corrected for the effect of the 
beak angle used in flying a turn by Interpolating the 
controls at an aircraft weight multiplied by the load 
factor 1/cos Integration step else varies during 
synthesis. During decelerations or accelerations It 
Is 1 sec while during altitude changes at fixed speed 
it is 5 sec. Total tima for synthesising s complete 
trajectory consisting of s horizontal trajectory 
similar to the ones shown in Fig. 3 end e speed/ 
altitude profile similar to tha ona in Fig. 5 la 
about 2 sac on the Sperry Flight Systems 1819A air- 


borne computer used in the flight tests. When the 
trajectory synthesis is rime-shared with navigation 
and either necessary computations, the computing time 
increases to about 6 sec. 

Real-Time Profile Generation 

After a profile has been synthesized in fast 
time and the pilot has elected to fly it, the refer- 
ence states and controls for that profile must be 
generated in real time. The synthesized profile 
can contain discontinuous changes in roll and pitch 
angle and throttle and vectoring nozzle position at 
’’command points” where changes in speed, altitude, 
or heading are initiated or terminated. The real 
time profile generation therefore must provide a 
certain amount of lead time to these control 
variables to minimize tracking errors at command 
points. These functions are performed by the Real- 
Time Profile Generation Logic. 

If an on-board computer has f ficient memory 
capacity to store all of the rete*. . e states and 
controls during fast ti*ne synthesis at small inter- 
vals of time, this logic would be relatively simple. 
However, limitations on the storage available in the 
STOLAND computer mode this approach impractical. To 
minimize memory usage, a different method was imple- 
mented at the expense of Increased complexity of 
computation. The method consists of storing refer- 
ence trajectory data, i.e., control positions, speed, 
altitude, etc., only at the "command points,” as 
defined earlier. Between "command points” the 
reference trajectory is generated in real time by the 
same integration logic used during fast time synthe- 
sis; however, the integration is now done entirely 
in forward time. Generation of a flyable reference 
trajectory that meets the desired boundary conditions 
is guaranteed because it is a precise repetition of 
a previously successful synthesis. 

The real-time forward integration uses disuiuue 
along the ground track as the independent variable. 
The integrated or dependent variables are refer- 
ence time, airspeed, altitude, and heading. The 
choice of distance rather than time gives a move 
flexible and operationally improved system for the 
following reasons. In a distance-based reference 
trajectory system, the aircraft will track the 
reference airspeed and altitude regardless of winds 
as it flies along the ground track. It is not 
necessary to null time errors if time control is not 
required. The system can thus be operated either in 
a 3D- or a 4D-tracking mode, depending on whether 
the time error loop is open or closed. This flexi- 
bility Is lacking in a time-based reference trajec- 
tory system, where only the 4D tracking mode is pos- 
sible. In the time-based system, if the actual 
winds differ significantly from the forecast winds 
used in fast-time trajectory synthesis, the air- 
craft controls may have insufficient authority to 
track the reference position resulting in unaccept- 
able tracking characteristics. 

One difficulty with the distance-based reference 
trajectory is that distance along the trajectory does 
not necessarily increase monotonicelly with time. 
Large navigation errors can cause tne new reference 
position to fell behind the previous one or to move 
ahead with e large step. This can result in control 
system saturation during the critical descent and 
deceleration segments. The system there 'ore contains 
logic that prevents the reference position from back- 
ing up or from advancing f after then about 1.5 times 
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the current ground speed between position updates, throttle, tt Cf falls in the engine operating range, 
which occur at 100 msec intervals. 84?. < tt < 9b%. 


Perturbation Guidance Law 


Perturbations of the aircraft states from the 
reference states are used in the guidance law to 
generate perturbation controls which are added to 
the reference controls in order to null errors in 
airspeed, altitude, and crosstrack position. The 
feedback states in the guidance lsw also include 
crosstrack error rate and flight-path angle as well 
as the integrals of airspeed and altitude errors. 

The latter two are used to reduce speed and altitude 
bias errors caused by inaccuracies in the stored 
energy-rate data and errors in the estimates of wind 
v id temperature profiles. 


The perturbation equations and the perturbat . :>n 
control law can be written in stat* vector notation 
as 

4r « Fx + Cu ; 14) 

dt 

u “ Kx (15) 

where 

x - (tv, £j, .\h, J'Vdt., /Ghdt ) 1 

, , . ,T 

u * ( ", ^v) 


The controls are throttle, nozzle, pitch, and 
roll angles. Flaps are not :$ed as perturbation 
controls because of their relatively low rate limit 
and an operational constraint that flap motion be 
monotonic during an approach. The flap command is 
simply the reference value at each ground track 
position limited to the placard value at the current 
airspeed. 


Lateral perturbation control is essentially 
uncoupled from the longitudinal mode and is accom- 
plished through a roll-angle command to the roll- 
command autopilot. This command is of the form 


f +k :y 
r %y 7 


+ 


where # $ r is the reference roll angle, and Ay 
and Ay are the crosstrack error and error rate, 
respectively. The two gains were chosen to provide 
a well-damped response and control activity compat- 
ible with the noise characteristics of the naviga- 
tion system. 


Longitudinal perturbation control for correct- 
ing airspeed and altitude errors is difficult 
because the refeiencc controls generated by the 
energy-rate schedule of Fig. 1 often lie on a con- 
straint boundary and therefore cannot be perturbed 
freely in both directions. The two controls that 
are often constraint limited during a fuel- 
conservative approach are throttle, and nozzle 
angle, v. Some Insight into this problem can be 
obtained using data from the energy-rate schedules. 
Figure 6 shows the energy-rate envelope from Fig. 1 
with the minimum reference nozzle and minimum refer- 
ence throttle constraint boundaries. These bound- 
aries divide the envelope into four regions: 

I, where v cannot be reduced; II, where neither 
it nor v can be reduced; III, where tt cannot be 
reduced; and IV, where it and v are free to move 
In either direction. The combinations Decontrols 
available for increasing and decreasing in 
each region are indicated in the figure. Note that 
in region I nozzle could be used as an additional 
control variable for decreasing energy rate. 

However, this variable is not used because throttle 
and pitch provide adequate control of flight-path 
errors in this region. In region IV the minimum 
reference throttle is above idle and Is determined 
by the maneuver margin constraint* At each air- 
speed in this region the negative throttle perturba- 
tion that can be added to the reference throttle to 
yield the commanded throttle is limited to -2% for 
safety reasons. Positive and negative throttle per- 
turbations art further limited so that the commanded 


The delta quantities are the perturbations from refer- 
ence values, i.e., AV - Y a - V ar , etc., where V a 
is the aircraft and Var the reference true air- 
speed, respectively. The commanded controls are 
the sum of reference and perturbation controls. 

u - (tt + At, e + AS v + A;) (16) 

c r 1 r 

For a powered-lift STOL aircraft such as the one 
used for these flight tests, the values of F and C 
are strongly dependent on airspeed and energy rate 
and are thus time-varying alcng a trajectory. 

Quadratic Optimal Synthesis" would therefore yield 
time-varying gain matrices that are also functions 
of the reference trajectory. But it is neither 
practical nor necessary to implement a complex, 
reference-trajectory-dependent gain matrix in order 
to achieve adequate control system performance ir* 
this case. 


The design procedure employed here began by 
first computing optimum gain matrices at various 
operating points in the control region diagram 
(Fig. 1) using fixed values of F and G. The 
analysis of these gain matrices showed the strongest 
dependence on airspeed, reference nozzle angle, and 
reference flaps. Sensitivity of the closed-loop 
eigenvalues to changes in several of the gains vac 
low, allowing those to be set to zero or held con- 
stant throughout the operating region. It was pos- 
sible to fit the variable gains with relatively 
simple functions of reference airspeed, nozzle angle, 
and flap angle. This method resulted In the fo 1 low- 
ing gain matrix: 

K - 



where V 4r is in unite of ft /see. Extensive 
computer calculations have verified that the closed- 
loop eigenvalues of this system have damping factors 
of 0.707 or greater and real parta less than -0.05/sec 
at all operating points. These characteristics 
provide adequate tracking performance. When 
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operating in region X of Fig. 6 the last row of K 
is set to zero since nozzle angle is not used for 
control. In regions II and III throttle pertur- 
bations are limited to positive values, while in 
region II nozzle perturbations are limited to posi- 
tive values. In region IV each control moves 
freely, but negative throttle perturbations are 
limited to -2% RPM, as previously explained. 

Control limiting can reduce the effectiveness of 
integral feedback of speed and altitude. Some 
design considerations for these integral feedback 
loops are given in Ref. 9. 

The throttle and nozzle angle perturbations 
generated by the control law will generally be of 
opposite sign, because the elements of the first 
row of K all have opposite sign of the third row 
elements. Thus, even in region II, where throttle 
and nozzle perturbations are each limited to move 
only in the positive direction; they are not gener- 
ally limited simultaneously. This implies that two 
controls, either throttle and pitch or nozzle and 
pitch, are free to move. Transient response studies 
using a nonlinear simulation of the aircraft and 
guidance system have shown that the control power 
fa adequate to provide rapid and well-damped air- 
speed and altitude error responses in region II. 
Example transient responses from this simulation 
are discussed in the last section. 

Structure and Operation of the Flight System 

The integration of the functional units of the 
system is shown in Fig. 7. Computations begin in 
the fast time trajectory synthesis module. If a 
trajectory Is successfully synthesized, It is stored 
at comand points, as previously explained, and the 
synthesized horizontal trajectory is display?! to 
the pilot on an electronic Horizontal Map Display 
(HMD) . This display operates, in conjunction wit 
the navigation system to give a map-like view of 
the terminal area (see Ref. 2 for details on this 
device). If the track switch is engaged, real-time 
profile generation and closed-loop tracking of the 
synthesized trajectory begins. The four control 
variables generated by the perturbation feedback 
iaw drive the roll and pitch autopilots and the 
throttle and vectoring nozzle servos* 

Figure 8 gives an example of trajectories dis- 
played on the HMD. The solidly drawn track is a 
fixed and prestored reference trajectory on which 
waypoint numbers are Indicated. The pilot selects 
the waypoint on the fixed reference trajectory he 
wishes ..o capture by keyboard entry (waypoint 2 in 
this case). The track drawn with broken lines from 
Pj to 2 Indicates to the pilot that a valid capture 
trajectory to that vaypofnt has been computed. If 
the synthesis had not been successful, the synthesis 
routine would have been reentered, as shown in 
Fig. 7, with updated aircraft states as the new 
initial conditions. 

To account for the distance the aircraft will 
travel while the trajectory is being synthesised and 
to give the pilot time to puah the track switch, ths 
capture trajectory is actually computed from Pj 
rather than from tha aircraft position at Pi (see 
Fig. 8). Tha short, straight segment between Pi 
and P 2 is rafarrad to as s lead distance and It ia 
drawn along tha aircraft velocity vector at tha 
start of synthesis. This technique minimises 
initial condition errors at tha beginning of tha 
h J tomatic-t racking mode. 


After the reference trajectory has been stored 
and valid navigation data from TACAN+ or KODILS* are 
received the track switch (see Fig. 7) is armed, 
ready for the pilot to engage. If the pilot does 
not engage che track switch before the aircraft 
reaches P 2 , P 2 i fi moved to its new position, and a 
new capture trajectory is displayed if one were suc- 
cessfully computed. When the pilot engages the 
track switch as at Pi in Fig. 8, the capture path 
on the display is drawn with a solid line indicating 
to the pilot that closed-loop guidance to the syn- 
thesized trajectory has begun. 

While the capture-tralectory algorithm synthe- 
sizes successful trajectories for a wide range of 
initial conditions, there are conditions where it 
will fail to do so. For example, if Pj is very 
close to the capture waypoint, then the algorithm 
can fail because there is insufficient distance 
along the computed minimum distance path to complete 
the required change In the speed and/or altitude. 

In that case, the reason for the failure to syn- 
thesize is displayed as a short message on the HMD. 
The pilot can correct the failure to capture condi- 
tion by flying the aircraft away from the capture 
waypoint or by selecting a more distant capture 
waypoint. 

The fixed reference trajectory, though not 
always usable as explained in the introduction, 
prescribes a nominal approach route and is determined 
by air-traffic control and terminal-area constraints. 
The precise airspeed and altitude profiles along 
this fixed reference usually are not rigldlv speci- 
fied. Often the airspeeds and altitudes are speci- 
fied only at waypoints. In that case the speed 
and altitude profiles between adjacent wavpoints 
are synthesized in fast time using the same algo- 
rithm as for the capture trajectory. The synthesis 
is done in backward time starting at the last way- 
point and ending at the capture waypoint. The 
altitude and speed at wavpolnt ?!-l determine the 
final condition and the altitude and speed at 
wavpolnt N determine the initial condition for 
the synthesis. Thus, every available degree of 
freedom is exploited to optimize the total 
t r<* Jec tory . 


Simulation and Flight Test Results 

The performance of the guidance system was eval- 
uated in simulation and flight tests. The piloted 
sinulator was the primary tool for determining the 
performance limits of the system since it allowed 
the measurement of performance for known disturbance 
inputs. In flight, it la difficult to measure or 
control disturbances and Isolate their effect on 
performance. Flight tests were used to verify the 
simulator model and to obtain pilot comments on the 
operational acceptability of the system. 

Figure 9(a) gives simulator time histories 
of selected stetes and controls for a combined cap- 
ture and fixed approach trajectory , from 3000 ft 
altitude, 140 KEAS, and 40,000 ft to touchdown. Tha 


i’TACtical Air Navigation System provides azimuth 
and range relative to station location. 

tHODu lar Instrument Lending System is an 
interim microwave landing system with a azimuth scan 
of t20*, an elevation scan of 16* , and a precision 
distance measuring system with a range of 10 n. mi. 
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Initial position and heading were chosen to yield 
a more cr less straight-in horizontal capture path 
along the extended centerline of the runway. Way- 
point 2 of the fixed-reference path shown in Fig. 8 
was selected as the capture waypoint. The capture 
trajectory consibts of a constant speed , level 
flight segment to point A where a -7.5* descent 
begins. Deceleration at 0.05 g starts at point B 
while the aircraft is in descent. The deceleration 
is initiated by starting the deployment of flaps 
(not shown in the figure). End of capture and the 
start of the f ina) -approach trajectory is at point 
C, which corresponds to waypoint 2 in Fig. 8, where 
the aircraft has decelerated to the landing speed ot 
72 KEAS and is tracking a -7.5° glide slope. 

Point C is 700 ft above the runway and about 5400 ft 
from touchdown. Automatic tracking is terminated at 
point 0 (Waypoint 1 in Fig. 8), 300 ft above the 
runway and 2280 ft from the nominal touchdown point. 
The capture trajectory was synthesized using a 
minimum flight-path angle of -7.5° and a maximum 
deceleration of 0.05 g. The parameters o and c 
were set to 0.9 and 1.0, respectively. The choice 
of e * 1 gives the priority to deceleration during 
the combined deceleration and descent segments. 

This combination of limits and parameter values 
resulted In a reduction in flight-path angle from 
-7.5* to -6.9* for a 10 sec period just prior to 
point C. Both the synthesis algorithm and the 
aircraft simulator assumed zero wind speed, a 
standard temperature profile, and an aircraft weight 
of 48,000 lb. Navigation errors were set to zero. 
The tracking performance of the perturbation- 
control law under these nominal conditions serves 
as a standard against which the performance under 
off-nominal and flight-teat conditions can ba 
evaluated. 

The step throttle reduction to 842 KPM and the 
pitch* down command (not shown in Fig. 9a) lead the 
computed descent point at point A by about 6 sac to 
compensate for throttle- and pitch-angl** dynamics. 
This prevents overshoots in tracking the reference 
altitude at the point of descent. The small alti- 
tude transient at this point is damped in about 
5 sec. Speed and altitude errors during the 
deceleration and descent segment between points B 
and C never exceed 4 knots and 20 ft, respectively 
These errors converge to zero between points C 
and D. Though of no practical significance, the 
residual errors during deceleration are caused by 
inaccuracies in the reference throttle nozzle, and 
pitch commands computed from the energy rate per- 
formance tables during fast time synthesis. The 
response of the perturbation control law to these 
modeling inaccuracies is seen as the difference 
between the superimposed traces of commanded and 
reference nozzle end throttle angle* in Fig. 9a* 
These differences are small even while reference 
nozzle and throttle angles are increasing rapidly 
toward tha end of the declaration segment. The 
validity of tha tnargy-rate-parformance model aa an 
accurate predictor of aircraft performance during 
these quasi-dynamic maneuvers is thus verified. 

Figure 9b shows tha r c..trol system response to 
a 10-knot headwind pulse of about 70-sec duration 
using the seme reference trajectory as in Fig. 9a. 
Tha initial spaed error resulting from this pulse 
is rapidly damped by using the nozsle ss control* 
Then, integral speed end altitude feedback develops 
throttle-, nozzle-, and pitch-parturbation biases in 
about 20 sec to cor r act the reference controls for 
tha change in aerodynamic flight-path angle caused 


by the change in the mean wind. In effect, the 
Integral feedback serves aa an estimator of the mean 
wind. When the headwind pulse is removed just after 
point C, the effect Is equivalent to a tailwind 
pulse and results in another transient. Note that 
the introduction of the headwind puJse hay increased 
the time between points B and C because the ground 
speed is reduced. CrosstracK position errors, 
though not shown, remain negligibly small during the 
approach. 

Other simulation results have shown that the 
control lav can compensate fc^ u nmodelid tailwinds 
and headwinds of 15 and 25 knots, respectively, with- 
out excessive control saturation. Also, errors of 
102 in aircraft weight and +20* C and -10° C 
unmode led temperature deviations from that assumed 
are compensated by the control lav. 

Figure 10 shows various time histories for a 
straight-in flight-test approach under conditions 
similar to those in the simulations. One significant 
difference was the use of .03 g maximum deceleration 
in the flight test. This cauced the deceleration to 
begin at point A before the point of descent at B. 

The entire approach trajectory was flown using TACAN 
for navigation. The location of the TACAN station 
relative to the runway Is shown In Fig. 8. Because 
of its favorable location, the TACAN station provided 
sufficient navigation accuracy for flying the 
approach without switching to the higher precision 
MODILS as would normally be required. There was 
light to moderate turbulence and an average head- 
wind of about 15 knots below 4000 ft as measured by 
a radar-tracked weather balloon just prior to take- 
off. However, the wind profile was not entered into 
the synthesis logic and thus constituted an unmodeled 
wind. Altitude errors, except near the pitch-down 
point, did not exceed 35 ft and decreased to about 
15 ft near the end. Speed errors during deceleration 
were less than 10 ft/tec, and decreased to about 
1 ft/aec at the end. If allowance is made for the 
presence of turbulence, winds, and navigation- system 
noise during the flight, these errors agree reason- 
ably well with the simulation results and are accept- 
ably low. The control perturbation biases, evidently 
caused by modeling errors and the unmodeled wind, are 
larger than those seen in the simulation run of 
Fig. 9a, though they are not excessive. Nozzle bias 
during the middle cf the deceleration averages about 
25*. While this seems large it should be noted that 
during this lntarval the throttle is at flight idle, 
where the effect of nozzle on energy rate is a mini- 
mum. On the whole, the control biases represent 
fairly small errors in the energy rate model. The 
flight- test results can, of course , be used to 
Improve the accuracy of the anargy-rate model of the 
aircraft. 

The crosetrack error at point D measured by 
precision redar wee 80 ft or less than half tha 
width of the runway. This error, which the pilot 
can null during manual 'light from point D (way- 
point 1) to touchdown, approaches tha accuracy 
limits of TACAN. At the same time it tats a lower 
limit on the distance between point D and touchdown 
whan using TACAN. Tilota judged the action of the 
control law as smooth and the capture trajectory ae 
a convenient end effective tool for optimizing 
approach trajectories. 

Tha fuel consumption of this automatically 
flown trajectory wee compared with that of a 
trajectory flown s y a test pilot under simulated 
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inatrument-f light-rule conditions. In order ti 
provide « bee is for comparison the manually flown 
trajectory began from the same Initial distance- tc- 
touchdown, airspeed and altitude as the automati- 
cally flown trajectory. The approach was made with 
the aid of a flight-director system which displayed 
to the pilot lateral and longitudinal deviations 
from a straight-in, 7.5* approach path. The fuel 
used for the automatic approach was 381 lb, while 
that for the manually flown one was 500 lb. Further 
simulations and flight tests should be conducted to 
compare the fuel consumption for various approach 
trajectories, fHght-dlrector designs, and wind 
conditions before the fuel conservation potential 
of this guidance system can be considered 
established. 


Concluding Remarks 

The automatic guidance system described in this 
paper achieves the dual goal of fully automatic 
flight and near-optimal fuel conservation through 
the technique of fast-time onboard trajectory syn- 
thesis. This technique overcomes the performance 
limitations Inherent in a stored, precalculated 
trajectory by adapting the trajectory to the unique 
conditions encountered in each landing approach. 

The ability to adapt is crucial in the terminal area 
since the initial conditions for starting the 
approach and the wind and temperature profiles are 
not predictable with sufficient accuracy prior to 
takeoff. The technique for synthesising the tra- 
jectories allows the pilot to choose various 
deceleration/descent profiles. All profiler have 
the comsun characteristic oi delaying the start of 
the descent and deceleration points as much as pos- 
sible. A preliminary flight evaluation indicates 
that an automatically flown, optimum approach can 
produce significant fuel savings relative to an 
approach flown manually with only conventional 
flight-director guidance. The design procedure 
described herein for s STOL slrcreft is applicable 
with lesser control complexity to guidance system 
design for CTOL elrcrsft. 
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Fig. 6 Constraint boundaries for perturbation controls. 
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Fig. 7 Block diagram of guidance system. 
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Fig. 8 Horizontal flight paths displayed on cockpit Horizontal Hap Display (HMD). 
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(b) Effect of 10-knot headwind impulse. 


Fig. 9 Simulation of typical approach trajectory; runway at 140-ft A.S.L 
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Fig. 10 Flight-test results, runway at lAO-ft A.S.L. 
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